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ABSTRACT 

Recent numerical relativistic results demonstrate that the merger of comparable-mass spinning 
black holes has a maximum “recoil kick” of up to ~ 4000kms _1 . However the scaling of these recoil 
velocities with mass ratio is poorly understood. We present new runs showing that the maximum 
possible kick parallel to the orbital axis does not scale as ~ rj 2 (where rj is the symmetric mass ratio) , 
as previously proposed, but is more consistent with ~ 77 s . We discuss the effect of this dependence on 
galactic ejection scenarios and retention of intermediate-mass black holes in globular clusters. 

Subject headings: black hole physics - galaxies: nuclei - gravitational waves — relativity 


1. INTRODUCTION 

Recently, numerical exploration of the radiative recoil 
“kick” of merging black holes has progressed consider- 
ably. In particular, efforts in this regard have led to sug- 
gested phenomenological formulae for the kick, largely 
based on post-Newtonian (PN) predictions such as that 
given by Kidder (1995), which have proved surprisingly 
successful. For example, Gonzalez et al. (2007) found 
that in cases of unequal masses ( q = mi/ m 2 < 1) and 
no spin, a simple modification of the PN formula orig- 
inally found by Fitchett (1983) fits the numerical data 
quite well. For cases of spins parallel or antiparallel to 
the orbital angular momentum, a formula proposed by 
Baker et al. (2007), loosely based on PN calculations, 
is also consistent with numerical data. For spins with 
components perpendicular to the orbital angular momen- 
tum, a formula, again derived from PN calculations, has 
been proposed by Campanelli et al. (2007). This formula 
agrees well with numerical results for equal masses. 

This last type of kick, which is parallel to the orbital 
angular momentum, is of particular interest because its 
computed magnitude can be very large (up to thousands 
of kilometers per second). However, unlike for the kicks 
that arise from unequal masses with no spins, or the in- 
plane kicks that arise from spins out of the orbital plane, 
the mass ratio dependence of these critical out-of-plane 
kicks has not been tested systematically by numerical 
experiments. In the current literature (specifically Cam- 
panelli et al. 2007), the dependence is drawn from the 
leading-order PN approximation. It is unclear whether 
this approximation is sufficient to predict the strong-field 
dynamics that presumably determines the kick. Indeed, 
hints of a deviation from this form are evident for mass 
ratio q = 1/2 in the runs of Lousto & Zlochower (2007). 
Therefore, although it can be shown that the angular de- 
pendence of the proposed formula is consistent with sym- 
metry arguments, which are independent of the strong- 
field dynamics (Boyle et al. 2007; Boyle & Kesden 2007), 
the mass ratio dependence of this formula is currently 
not well justified. 
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Characterization of the dominant kick for unequal 
masses is especially important because, although the 
largest possible kicks would eject the remnant from any 
galaxy, for astrophysical applications it is the distribu- 
tion of kick speeds that matters. For example, Bonning 
et al. (2007) find no evidence for quasars ejected from 
their hosts. If quasar activity is commonly induced by 
major galaxy mergers that lead to coalescence of super- 
massive black holes, the implications of this therefore 
depend in part on how frequently one expects a merger 
to allow ejection. In addition, the kick speed distribution 
has a major impact on the hierarchical growth of massive 
black holes at redshifts z > 5 (e.g., Volonteri 2007). 

Here we investigate how the out-of-plane kick depends 
on the mass ratio, and find that, for mass ratios in the 
range <7 = 1 to <7 = 1/3, the kick drops off more rapidly 
with decreasing mass ratio than proposed by Campan- 
elli et al. (2007). Specifically, we find that all current 
numerical data on kicks are well represented by 

Vrecoi\=Vmh + vj.(cos£ei + sin £ e 2 ) + ^|| e 2 , (1) 

v m = ATi 2 y/i -47?(1 + B?7), (2) 

(3) 

Krj 3 

W|| = ^ [qai cos (4>i - $i) - ot% cos(<p 2 - $2^4) 

where r\ = qj (1 + q) 2 is the symmetric mass ratio, af 
is the projection of the dimensionless spin vector Si = 
Si / m ? of black hole i along the orbital angular momen- 
tum, af- is the magnitude of its projection, af~, into the 
orbital plane, fa refers to the angle made by a f~ with re- 
spect to some reference angle in the orbital plane, and 4>i 
and $2 are constants for a given mass ratio. Here, A = 
1.35 x 10 4 kms -1 , B = -1.48, H = 7540 ± 160 km s~\ 
£ - 215° ±5°, and K = 2.4 ± 0.4 x 10 5 kms" 1 . This for- 
mula, similar in form to that given by Campanelli et al. 
(2007), synthesizes results from Gonzalez et al. (2007) for 
(2) and from Baker et al. (2007) for (1) and (3) 3 . For 
£ and H we have fit available numerical data from Her- 
rmann et al. (2007); Koppitz et al. (2007); Baker et al. 

3 Note that in Baker et al. (2007), we used a simpler form for 
the zero-spin contribution, equivalent to (2) with B = 0. 
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TABLE 1 

Initial Configuration and Final Kick for Each 
Simulation. <p i( 2 > is the angle made by the spin 

VECTOR OF HOLE 1(2) WITH THE VELOCITY VECTOR OF 
HOLE 1, AS SHOWN IN FlG. 1. NUMERICAL RESULTS FOR 
THE KICK COMPONENTS V m (WHERE AVAILABLE) AND 1>|| 

are shown. Kicks for equivalent spinless runs 

ARE IN PARENTHESES. 


q 

M° ) 

M°) 

v m (kms 1 ) 


( km s 1 ) 

1/1.1 

0 

180 

24 


-542 


315 

135 

24 


-657 


270 

90 

25 


-384 

1/1.3 

0 

180 

67 


-386 


315 

135 

67 


-525 


270 

90 

69 


-348 

1/1.5 

60 

240 

92 (94) 


-381 


0 

180 

95 (94) 


-135 


315 

135 

91 (94) 


168 


270 

90 

90 (94) 


364 

1/2 

0 

180 

137 (140) 


-37 


315 

135 

136 (140) 


111 


270 

90 

136 (140) 


193 


315 

90 



75 


0 

90 



-55 

1/3 

0 

180 

166 


49 


315 

135 

166 


48 


270 

90 

163 


17 

(zuu/j. ine qualitatively new part, 

the 

factor 6f 


(4), replaces the factor of rj 2 originally proposed by Cam- 
panelli et al. (2007), and is motivated by new numerical 
evolutions presented here. We give our methodology in 
§ 2, and present our results, fits and possible analytical 
motivation, and astrophysical consequences in § 3. 

2. INITIAL DATA AND METHODOLOGY 

We simulated the inspiral and merger of a range of 
spinning black-hole binaries, with mass ratios in the 
range 1/1.1 > q > 1/3. The initial configuration of 
momenta and spins is illustrated in Fig. 1. The parame- 
ters used in the numerical evolutions are presented in the 
first three columns of Table 1. For these evolutions, the 
smaller hole (mi) has a dimensionless spin |<5i| = 0.2, 
while the larger hole’s spin is |c?2| = Q 2 |cfi|- Both spin 
vectors initially lie in the orbital plane, at angles (pi and 
02 to the initial velocity of hole 1 (see Fig. 1). 

To perform our simulations, we employed the Hahn- 
DOL evolution code, as described in Baker et al. (2007). 
The code’s convergence properties are discussed in Baker 
et al. (2007) and in Baker et al. (2008). Initial separa- 
tions were chosen to yield between one and four orbits 
prior to merger; the corresponding momenta were cho- 
sen, informed by PN theory (Damour et al. 2000), to 
minimize initial eccentricity. 


3. RESULTS AND DISCUSSION 

The recoil kicks resulting from the new simulations are 
given in the rightmost columns of Table 1. Now we con- 
sider the agreement of our data with the suggested for- 
mula of Campanelli et al. (2007): 

r? 

V|| = Kcos(@ - @o)jY—j\®2 -qail ( 5 ) 



Fig. 1. — Configuration of black holes for all new simulations. 
The two holes’ spins <3( 12 ) lie initially in the orbital plane, at angles 
(pi and <p2 to v\ , the smaller hole’s initial velocity. 

where 0 is the angle between (qa i — a^) and the sep- 
aration vector, nominally just before merger, and 0 o is 
some constant for a given mass ratio. In practice, it is 
more convenient to consider © with respect to the ini- 
tial separation vector, and absorb into ©o the amount 
by which this angle precesses before merger. By this def- 
inition, then, ©o depends on both the mass ratio and 
the initial separation. Also note there is a second angle 
implicit within the expression |a 2 — qai |, which repre- 
sents the angle between and a^- Recognizing that 
cos(© — ©o)|d?2 — qot i| = (a?2 — qa i) • n, where h is the 
initial unit separation vector, we can rewrite (5) as 

rj 2 

v\\ — K— - — r[qa± cos(0i-$i)-a^cos(0 2 -$2))] (6) 

(l + <7) 

where (pi represents the initial angle of af, and $1 and 
$2 are parameters which depend on mass ratio and initial 
separation, interpretable as encoding both the spin ori- 
entations for the maximum kick as well as the amount of 
spin precession before merger. Here we keep the value of 
K — 6.0 x 10 3 4 , which was found to work well in the equal- 
mass case, and find the best fit for the parameters $1 and 
$2 for each mass ratio. The error of the best fit grows sig- 
nificantly with mass ratio (as seen in the column labeled 
by “Krj 2 ” in Table 2), hence the mass-ratio-dependence 
of this formula is inaccurate. One might suppose that 
precession of the spins into the orbital plane could ac- 
count for this. However, the v m column in Table 1 shows 
that the in-plane kicks are close to those measured with- 
out spins (given in parentheses), hence this is not an 
explanation of the discrepancy in U|j from the r ] 2 scaling. 

To conceive of other plausible candidates for the kick 
formula we begin with the spin expansion and symmetry 
arguments of Boyle et al. (2007); Boyle Kesden (2007). 
For the spin configurations under consideration here, 

i’ll = D(q)a ^ cos (0i - $(q))-D(l/q)oi 2 cos(0 2 - $(l/g)) 

(7) 

where D and are some functions of mass ratio #, and we 
note that $ must also depend on the initial separation. 
Further restricting ourselves to forms relatable to the 
factor of S\ /mi — S^/m^ appearing in PN calculations 
of the kick, which informed Campanelli et al. (2007); 
Lousto & Zlochower (2007) and has been numerically 
well- verified in the equal-mass case, we substitute D(q) = 
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TABLE 2 

Maximum percent error resulting from various models 

OF THE KICK, AS DISTINGUISHED BY OVERALL MASS-RATIO 
DEPENDENCE. SEE EQUATION (8). 


Q 1 

Krj 2 

K(a2T] 2 4 * CI477 4 ) 

Kr ) 3 

1/1.1 

0.22 

0.23 

0.20 

1/1.3 

0.75 

0.80 

0.78 

1/1.5 

1.26 

1.31 

1.28 

1/2 

15.57 

2.46 

1.41 

1/3 

39.37 

10.81 

1.82 


qC(rj)/( 1 +q) to obtain: 

«|| = [l 01 ! COS (<l>l - $l) - a 2 C0S(^ 2 - $2)] , 

q (8) 

where 4>i = $(q) and $2 = 4>(l/g). 

Regarding the overall scaling of v\\, it is known to 
be related to the difference between the energy radi- 
ated in the (Z,ra) = (2,2) and (2,-2) harmonics of 
the radiation (Briigmann et al. 2007). With no spin, 
these quantities are equal. With spin, we expect that 
V\\ ~ E 2 2{peak)F , where E 22 ( P eak) is the peak energy ra- 
diated in the (2, 2) harmonic, and F represents the spin- 
dependent asymmetry between E 2 2 ( peak ) and E 2 _ 2(<peak) , 
i.e. F ~ 1 - E 2 __ 2{peak) /E 22 ( peak ). For black holes with 
no spin it has been found that E 22 ^ peak ) = a 2 q 2 -f a^r] 4 , 
where a 2 = 0.0044 and a\ = 0.0543, gives a good fit to 
the numerical data. We do not expect spins orthogonal 
to the orbital angular momentum to change the scaling 
of the radiated energy significantly. If we further assume 
that the asymmetry factor F is independent of 77, which 
finds some support in PN analysis since to leading order 
P\\/E is independent of 77, then we are led to hypothesize 
that m oc (a 2 r] 2 + a^rj 4 ). 

In Table 2 we summarize the agreement of various kick 
formulas with the numerical data. For each formula, 
which has the form of Eq. (8), we found the best <3>i 
and <$2, per mass ratio, according to a least squares fit 
to the data given in Table 1. For each mass ratio, the re- 
sulting percent error is given for each model, maximized 
across initial angle. Referring to Eq. (8), the column 
headings Kr f, K(a 2 r f + a^q 4 ) and Krf* of Table 2 rep- 
resent choices for C(q) that were tested, where in each 
case K has been chosen so as to reproduce the value 
of the formula of Campanelli et al. (2007) in the equal- 
mass case. Note that C(rj) = Kr] 2 gives exactly the for- 
mula of Campanelli et al. (2007). We see that the choice 
C{rj) = K{a 2 r] 2 -b CZ477 4 ), for which we have strong-field 
heuristic justification, fits the data much more success- 
fully than does C(q) = Kr] 2 . However, a better empirical 
model was found to be C(rj) = Krf. For now we consider 
this our best fit, and leave open the interesting question 
of how to accurately relate this prefactor directly to E 22 . 

Our results affect the distribution of kick speeds given 
various assumptions about the spin parameters, spin ori- 
entations, and mass ratios involved in coalescences. This 
has particular application to the retention of the prod- 
ucts of mergers of massive black holes in the current uni- 
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verse (e.g., Bonning et al. 2007) and electromagnetic sig- 
natures of kicks (e.g., Shields et al. 2007; Lippai et al. 
2008), as well as coalescences in the early structure for- 
mation phase of redshift z ~ 5 — 30 (Merritt et al. 
2004; Boylan-Kolchin et al. 2004; Haiman 2004; Madau 
& Quataert 2004; Yoo & Miralda-Escude 2004; Volon- 
teri & Perna 2005; Libeskind et al. 2006; Micic et al. 
2006; Volonteri 2007), and for current-day mergers of 
intermediate-mass black holes (IMBHs), which might ex- 
ist in dense stellar clusters (Taniguchi et al. 2000; Miller 
& Hamilton 2002b, a; Mouri & Taniguchi 2002b, a; Miller 
& Colbert 2004; Giiltekin et al. 2004, 2006; O’Leary et al. 
2006, 2007). Note that q = 1 to q = 1/3, is in the range 
of ratios expected for major mergers of galaxies, and as 
Sesana et al. (2004) show, this range is expected to ac- 
count for most massive black hole mergers in the early 
z > 10 phase of black hole assembly. 

Our new formula implies an important revision in our 
understanding of how easily IMBHs with M ~ 10 2 — 
10 3 M© are retained in globular clusters. A rich cluster 
has an escape speed v esc « 50 kms -1 (Webbink 1985). 
Giiltekin et al. (2006) showed that the Newtonian kicks 
involved in binary-single interactions are insufficient to 
reach this speed if the IMBH is at least ~ 15 — 20 
times more massive than the objects with which it in- 
teracts. Using the Campanelli et al. (2007) formula, 
however, the maximum kick from gravitational radiation 
is u max = 6 x 10 4 kms -1 77 2 , implying that even IMBHs 
30 — 35 times more massive than the black holes with 
which they merge could get ejected. Holley-Bockelmann 
et al. (2007), focusing on cases in which stars lose little 
mass through their evolution and thus can leave behind 
stellar-mass black holes with masses > 60 — 100 M©, use 
this to argue that most IMBHs of even 1000 M© will 
be ejected from globulars. If instead stellar-mass black 
holes have masses ~ 10 M©, a mass of at least 400 M© 
would still be required to guarantee retention. 

In contrast, our new formula suggests a maximum kick 
of u max = 2.4 x 10 5 kms“ 1 77 3 . Thus if 77 < 0.06, u max < 
50 kms -1 . Therefore, an IMBH interacting with 10 M© 
black holes will stay in a rich globular if its initial mass 
is M > 170 M©, comparable to what is necessary for 
retention against Newtonian three-body kicks. 

Our results also have implications for whether merged 
supermassive black holes stay in their host galaxies. The 
figure of merit is the fraction of kicks that exceed typ- 
ical escape speeds from galactic centers (ranging from 
roughly 500 kms" 1 for a small spiral to 2000 kms -1 
for a giant elliptical), given assumptions about the dis- 
tribution of spins and orbital orientations. The calcula- 
tion of record for this purpose is that by Schnittman & 
Buonanno (2007), who used a kick formula based on ef- 
fective one-body analysis and is different from that of 
Campanelli et al. (2007); this formula underestimates 
the highest kicks. Table 3 compares the fraction of kicks 
above 500 kms -1 and 1000 kms -1 using the Schnittman 
& Buonanno (2007) formula (an underestimate), the 
Campanelli et al. (2007) formula (an overestimate), and 
our results. It is clear that the Schnittman Sz Buo- 
nanno (2007) results were conservative: the fraction of 
large kicks is significantly higher than their estimate for 
comparable-mass mergers with plausible spins. 

One consequence of the higher kicks is that retention 
of supermassive black holes after galactic major mergers 
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TABLE 3 

Fraction of kick speeds above a given threshold, compared with the results of Schnittman Sz 
Buonanno (2007) (SB) and Campanelli et al. (2007) (CLZM). In all cases we assume an isotropic 

DISTRIBUTION OF SPIN ORIENTATIONS. 


Mass ratio and spin 

Speed threshold 

SB 

CLZM 

This work 

1/10 < q < 1, a\ — a ,2 = 0.9 
1/4 < 9 < 1, ai = 02 = 0-9 
1/4 < q < 1, 0 < ai,a 2 < 1 

v > 500 km s -1 
v > 1000 km s -1 
v > 500 km s -1 
v > 1000 km s -1 
v > 500 km s -1 
v > 1000 km s -1 

0.12(+0.06, -0.05) 
0.027(+0.021, -0.014) 
0.31(+0.13, -0.12) 
0.079(+0.062, -0.042) 

0.364±0.0048 

0.127±0.0034 

0.699±0.0045 

0.364±0.0046 

0.428±0.0045 

0.142±0.0034 

0.2283±0.0014 

0.085±0.0008 

0.618±0.0014 

0.2547±0.0013 

0.3484±0.0015 

0.0974±0.0009 


is even more challenging than previously thought, un- 
less an astrophysical mechanism restricts the spin mag- 
nitudes (contrary to spin inferences from Fe Ka lines; 
see Iwasawa 1996; Fabian et al. 2002; Reynolds & Nowak 
2003; Brenneman & Reynolds 2006) or the spins tend 
to align parallel to each other and to the orbital axis 
(Bogdanovic et al. 2007). Absent such a mechanism, one 
would expect tens of percent of merged galaxies to have 
no central black hole, in strong contradiction with obser- 
vations (see Ferrarese & Ford 2005). Given that purely 
gravitational precession and radiation do not preferen- 
tially align spins in weak gravity (Schnittman 2004; Bog- 
danovic et al. 2007), nor are they expected to in strong 
gravity (A. Buonanno, private communication), align- 
ment would have to come from external torques, e.g., by 
nuclear gas if there is a sufficient amount in the vicinity. 

In conclusion, we have performed a systematic study 
of the mass ratio dependence of the out-of-plane kicks 
produced by the merger of spinning black holes. Our 
work shows that the Campanelli et al. (2007) candidate 


kick formula overestimates the out-of-plane kick system- 
atically. However, we find that an additional factor of 
Arj agrees with our numerical results to within typical 
values of 1% for mass ratios between 1 and 1/3. This 
has considerable implications for black hole retention in 
early dark matter halos, galaxies, and globular clusters. 
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